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Z9* Production in Deep Inelastic Scattering 

at HERA 


ZEUS Collaboration 


Abstract 

This paper presents measurements of D*± production in deep inelastic scattering 
from collisions between 27.5 GeV positrons and 820 GeV protons. The data have been 
taken with the ZEUS detector at HERA. The decay channel D* + —> (D° —> K~tt + )'k + 
(+ c.c.) has been used in the study. The e + p cross section for inclusive D * ± produc¬ 
tion with 5 < Q 2 < 100 GeV 2 and y < 0.7 is 5.3 ± 1.0±0.8nb in the kinematic region 
1.3 < < 9.0 GeV and | \ < 1.5. Differential cross sections as functions of 

P'Y^D*^), W and Q 2 are compared with next-to-leading order QCD calculations 

based on the photon-gluon fusion production mechanism. After an extrapolation of the 
cross section to the full kinematic region in p^^D^) and the charm contribution 

F£ c (x,Q 2 ) to the proton structure function is determined for Bjorken x between 2- 1CD 4 
and 5 • 1CT 3 . 
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1 Introduction 


In neutral current deep inelastic scattering (DIS), ep —> eX, charmed quarks are expected to be 
produced predominantly via the photon-gluon fusion (PGF) process which couples the virtual 
photon to a gluon of the proton. The leading order (LO) diagram is shown in Fig. |l]a. 

Recently, analytic calculations of the DIS charm cross section from photon-gluon coupling have 
become available [|l|, which relate the DIS charm cross section to the gluon distribution in the 
proton using next-to-leading order (NLO) QCD. 

Measurements of deep inelastic neutral current scattering at HERA have demonstrated a rapid 
rise of the proton structure function F 2 as Bjorken-x decreases below 1CT 1 2 J2J. A QCD analysis 
of these data has connected this rise to an increase of the gluon momentum density in the 
proton with a dependence x~ x with A = 0.35 Q 2 = 7 GeV 2 || f|. At HERA, the 

DIS charm cross section is sensitive to the gluon distribution of the proton at low fractional 
momentum (x g ~ 10~ 3 ) of the proton. Comparison of the DIS charm cross section with NLO 
QCD calculations allows an independent check of the increase of the gluon momentum density, 
testing the consistency of the QCD calculations. 


The other processes, apart from PGF, that contribute to open charm production in DIS are: 
diffractive heavy flavour production j5j, scattering of the virtual photon off the charm sea quark 
H , charmed hadron production from bb [|7| and production of cc in fragmentation || . These 
processes, however, are expected to have much smaller cross sections than PGF. The possible 
contribution from intrinsic charm |9j. || is outside the acceptance of the main detector. 


ZEUS [|1(| [UJ and HI |T2| have reported on D* production^ by quasi-real photons at HERA 
(i.e. in the photoproduction regime). In JlCj first D* signals in the DIS regime were also shown. 


In this paper a detailed study of the D* produced in DIS events is presented. D* are investigated 
in the decay channel 

D* + -> D° 7T+ -> iP"7r + 7r+(+ c.c.) (1) 

using a procedure first proposed in |H|, where 7T+ stands for the ‘soft pion’. 


The e + p cross section for inclusive D*± production and differential cross sections as functions 
of pi(D*), rj(D*), W and Q 2 are presented. The measurements are compared with LO and 
NLO QCD analytic calculations JIJ] based on the photon-gluon fusion production mechanism. 
By extrapolation, the charm contribution to the proton structure function F 2 , F| c (x, Q 2 ), is 
estimated and compared with the NLO QCD analytic calculations. A similar analysis has been 
recently presented by the HI Collaboration |b| . 


2 Experimental setup 


The data were collected at the positron-proton collider HERA using the ZEUS detector during 
the 1994 running period. HERA collided 27.5 GeV positrons with 820 GeV protons yielding a 

1 In this paper charge conjugate modes are always implied. D* always refers to both D* + and the charge 

conjugate mode (c.c.) D*~. 
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center-of-mass energy of 300 GeV. 153 bunches were filled for each beam, and, in addition, 15 
positron and 17 proton bunches were left unpaired for background studies. The r.m.s. of the 
vertex position distribution in the Z direction^ is 12 cm. The data used in this analysis come 
from an integrated luminosity of 2.95 pb _1 . 


ZEUS is a multipurpose detector which has been described in detail elsewhere [15]. The key 


component for this analysis is the central tracking detector (CTD) which operates in a magnetic 
field of 1.43 T provided by a thin superconducting solenoid. The CTD is a cylindrical drift 
chamber consisting of 72 layers organized into 9 “superlayers” covering the polar angular region 
15° < 6 < 164° m • Five of the superlayers have wires parallel (axial) to the beam axis and 
four have wires inclined at a small angle to give a stereo view. The spatial resolution in the 
drift direction is 190 /jrn. The interaction vertex is measured with a resolution of 0.4 cm in the 
Z direction and 0.1cm in the XY plane. The momentum resolution for tracks traversing all 9 
superlayers is ct(pt)/pt = 0.005 p^ © 0.016 (pr in GeV). 

The solenoid is surrounded by a high resolution uranium-scintillator calorimeter (CAL) de¬ 
scribed elsewhere [|T7j . The position of positrons scattered close to the positron beam direction 
is determined by a scintillator strip detector (SRTD) Jl 8 [] . The luminosity is measured via the 
Bethe-Heitler process, ep —► ep r 7 , where the photon is tagged using a lead-scintillator calorime¬ 
ter [1191 located at Z = — 107 m in the HERA tunnel. 


3 HERA kinematics 


The kinematics of deep inelastic scattering processes at HERA, e + [k) + p(P) —■> e + {k') + A", 
where X is the hadronic final state, can be described by the Lorentz invariant variables Q 2 , x and 
y. Here — Q 2 is the square of the four-momentum transfer at the lepton vertex, x is the Bjorken 
variable and y is the fractional energy transfer between the positron and the proton in the 
proton rest frame. In the absence of QED radiation, Q 2 = —q 2 = — (k — k') 2 , and x = 


Q 


2Pq ’ 


where k and P are the four-momenta of the incoming particles and k' is the four-momentum of 
the scattered positron. The variables are related by Q 2 = sxy , where s is the squared invariant 
mass of the ep system. Since the ZEUS detector is nearly hermetic, for neutral current DIS 
Q 2 , x and y can be calculated from the kinematic variables of the scattered positron, from the 
hadronic final state variables, or from a combination of both. In this paper we use the double 
angle method (DA) 0 to calculate the Q 2 , x, y variables. The center-of-mass energy of the 
virtual photon-proton system (y*p), W, is determined using W^ A = m 2 + Q 2 L >a(^J ~ 1), 
being the proton mass. The variable y, determined using the energy E' e and angle 6' e of the 
scattered positron, y e — 1 — (1 — cos#'), is used for background suppression. The use of 

the positron information alone to calculate Q 2 , x and y is used in the study of the systematic 
uncertainties. 


Using the calorimeter information the quantity 5 = J2i{E t — Pz t ) is measured, where E { is 
the energy and pz t the longitudinal momentum assigned to the calorimeter cell i. For perfect 
detector resolution and acceptance, S = 55 GeV for DIS events while for photoproduction 

2 The ZEUS coordinate system is defined as right-handed with the Z axis pointing in the proton beam 
direction and the X axis horizontal pointing towards the center of HERA. The polar angle 0 is defined with 
respect to the positive Z direction. 
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events and events with hard initial state radiation, where the scattered positron or the radiated 
photon escapes down the beam pipe, 5 peaks at lower values. 


4 DIS event selection 


4.1 Trigger selection 


The trigger selection is almost identical to that used for the measurement of the structure 
function F 2 [^1J. Events are filtered online by a three level trigger system [^2[. At the first 
level DIS events are selected by requiring a minimum energy deposition in the electromagnetic 
section of the CAL. The threshold depends on the position in the CAL and varies between 3.4 
and 4.8 GeV. At the second level trigger (SLT), beam-gas background is further reduced using 
the measured times of energy deposits and the summed energies from the calorimeter. The 
events are accepted if 5 calculated at the SLT level using the nominal vertex position satisfies 
$slt > 24 GeV — 2E 7 , where E 1 is the energy deposit measured in the luminosity photon 
calorimeter. 


The full event information is available at the third level trigger (TLT). Tighter timing cuts as 
well as algorithms to remove beam halo muons and cosmic muons are applied. The quantity 
5tlt is determined in the same manner as for Sslt ■ The events are required to have 5tlt > 
25 GeV —2 E 1 . Finally, events are accepted as DIS candidates if a scattered positron candidate 
of energy greater than 4 GeV is found. For events with the scattered positron detected in the 
calorimeter, the trigger acceptance is essentially independent of the DIS hadronic final state. 
It is greater than 90% for Q 2 « 5 GeV 2 (lower limit for the DIS sample in this analysis) and 
increases to 99 % for Q 2 > 10 GeV 2 as determined from MG simulation. 


4.2 Offline event selection 


The selection of DIS events is similar to that described in our earlier publication [^1]]. The 
characteristic signature of a DIS event is the scattered positron detected in the uranium scin¬ 
tillator calorimeter. The positron identification algorithm is based on a neural network using 
information from the CAL and is described elsewhere The efficiency for finding the scat¬ 
tered positron is sensitive to details of the shower evolution, in particular to energy loss in the 
material between the interaction point and the calorimeter. The efficiency of the identification 
algorithm when the scattered positron has an energy of 8 GeV is 50%, rising to 99 % for energies 
above 15 GeV. For the present data, 89 % of the events have an e + energy greater than 15 GeV. 
The impact position of the positrons is determined by either the position reconstructed by the 
CAL or by the SRTD if inside its fiducial volume. The resolution of the impact position is 
about 1cm for the CAL and 0.3cm for the SRTD. The positron traverses varying amounts of 
material in the detector before entering the CAL, which causes a variable energy loss. This en¬ 
ergy loss is corrected on an event by event basis, as explained in [|T8|. The scale uncertainty on 


the energy of the scattered positron after these corrections is 2% at 10 GeV linearly decreasing 
to 1% at 27.5 GeV. 
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The following criteria are used to select DIS events: the presence of a scattered positron candi¬ 
date with a corrected energy E' e > 8 GeV, and 5 > 35 GeV to remove photoproduction events 
and to suppress events with hard initial state radiation. After applying these criteria, 4.3 • 10 5 
events are retained. For this analysis the region 5< Q 2 <100 GeV 2 and y e < 0.7 is selected 
which contains 3.7 • 10 5 events. A subsample of DIS events in the (x, Q 2 ) plane is plotted in 
Fig. 0b. The lines of constant y and Q 2 delimiting the region chosen for this study are shown. 
The photoproduction contamination is less than 2% and the beam gas background is negligible. 


4.3 Event selection efficiency 


The efficiencies of the event selection and the D* reconstruction are determined using a GEANT 
based Monte Carlo (MC) simulation program which incorporates the knowledge of the 
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detector and the trigger. For this analysis two different types of event generators were used: 
neutral current (NC) DIS generators, discussed below, for event selection efficiency calculations 
and heavy flavour generators, discussed in section 5.2, for D* reconstruction efficiency estimates. 

Neutral current DIS events with Q 2 > 4 GeV 2 were generated using the HERACLES 4.4 


which incorporates first order electroweak corrections. 

m 


to either LEPTO 6.1-3 E7I1 or ARIADNE 4.03 


It was interfaced using 
for the simulation of 


program [[ 

DJANGO 6.1-2 

QCD cascades. The calculation of the zeroth and first order matrix elements plus the parton 
shower option (MEPS) was used in LEPTO. The latter includes coherence effects in the final 
state cascade via angular ordering of successive parton emissions. In ARIADNE, the colour- 
dipole model including the photon-gluon fusion process (CDMBGF) was used. In this model 
coherence effects were implicitly included in the formalism of the parton cascade. The Lund 
string fragmentation |29|], as implemented in JETSET 7.4 
phase. 


30 , was used for the hadronisation 


parameterisation was used for the MEPS data set. For the CDMBGF event 
[3^j parton density parameterisation for the proton was used. These 


The GRV-HO 
sample the MRSD—' 

parameterisations describe the HERA measurements of the proton structure function F 2 
p] reasonably well. 
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Monte Carlo samples of DIS events containing D* mesons (all decay modes) corresponding to 
5.5pb _1 were generated. The CDMBGF sample was used for cross section determinations and 
the MEPS sample for systematic studies. The efficiency for selecting DIS events where a D* 
has been produced is determined as a function of Q 2 and y. For the kinematic region selected, 
an average event selection efficiency of 75% is found, increasing from about 60 % to 90 % as Q 2 
increases from 5 to 100 GeV 2 and from about 65 % to 85 % as y increases from 0 to 0.7. 


5 D* Reconstruction 


The tracks of charged particles are reconstructed using the CTD. The single hit efficiency of 
the chamber is greater than 95%. The efficiency for assigning hits to tracks depends on several 
factors, for example the px of the track and the number of nearby charged particles. In addi¬ 
tion the 45° inclination of the drift cells, which compensates for the Lorentz angle, introduces 
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some asymmetry in the chamber response for positive and negative particles, particularly at 
low pr- The reconstructed tracks used are required to have more than 20 hits, a transverse 
momentum p t > 0.125 GeV and a polar angle between 20° < 9 < 160°. In terms of pseudora¬ 
pidity, rj = — log(tan(6 l /2)), this angle corresponds to \r)\ < 1.75. This is the region where the 
tracking detector response and systematics are sufficiently understood. For those tracks with 
Pt > 0.125 GeV and | ?/1 < 1.75 the track reconstruction efficiency is greater than 94%. 


5.1 D* Identification 

D* production is investigated in the decay channel (1). The tight kinematic constraint on the 
D* + —> D°7r+ decay limits the momentum of the decay products to just 40 MeV in the D* rest 
frame. This fact allows one to measure the mass difference M(D*) — M(D °) more accurately 
than the measurement of the D* mass itself. In practice it leads to a prominent signal in the 
AM = M^KiriVs) — M(Kn) distribution, in an otherwise highly suppressed region of phase 
space. 

The D* reconstruction procedure consists of two steps. First a D° pre-candidate is formed by 
taking all combinations of pairs of oppositely charged tracks and assuming each track in turn to 
correspond to a kaon or a pion. If the Kn invariant mass of the track combination lies between 
1.4 and 2.5 GeV, the track pair is considered to be a D° pre-candidate. In the second step this 
D° pre-candidate is combined with a third track, which has the sign of charge opposite to that 
of the kaon candidate. The third track is assumed to be a pion (the so called soft pion, n s ). If 
the mass difference, AM, is below 180 MeV, the three tracks form a D* pre-candidate. 

To reduce the combinatorial background and restrict the analysis to a kinematic region where 
the detector reconstruction efficiency is acceptable, the following requirements are applied: 


• The spatial resolution of the CTD does not allow the D° decay vertex to be distinguished 
from the primary vertex. Therefore, we require all tracks to be associated with the 
reconstructed primary vertex of the event. 

• The kaon and pion candidates from the D° decay must have transverse momenta greater 
than 0.4 GeV. This reduces the combinatorial background in the D° reconstruction step. 

• D* pre-candidates must have transverse momenta in the range 1.3 < p^ < 9.0 GeV and 
directions of flight away from the beams, |?/| < 1.5. The lower p x limit is due to the very 
small acceptance for D* with p x below 1.3 GeV which results from the px( 7r s)>0.125 GeV 
cut. The higher p ? limit is due to the lack of statistics. Note that these cuts limit the 
actual acceptance to y > 0.015. 

The analysis is restricted to those combinations (called D* candidates) which pass the above 
requirements, have a Kn invariant mass in the range 1.8 < M(Ktt) < 1.92 GeV and have a 
mass difference in the range 143 < AM < 148 MeV. Fig.[l]b shows the distribution of all DIS 
events which have a D* candidate in the (x, Q 2 ) plane. 

Figs.|2a,b show the resulting AM (M(K 7r)) spectrum for those pre-candidates with M(Ktt) 
(AM) inside the corresponding signal region. Clear signals above the combinatorial background 
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are seen around the nominal M(D*) — M(D°) and M(D°) values. The AM spectrum is fitted 
with the maximum likelihood method by a Gaussian shaped signal plus a background of the 
form dN/dAM = a (AM — m n ) b in the 140-180MeV mass region where a and b are free 
parameters. This fit yields M(D*) — M(D°) — 145.44 ± 0.09MeV, in good agreement with the 
PDG value of 145.42 ± 0.05 MeV. The width of the signal is 0.65 ± 0.10 MeV. The excess 
of events seen in Fig. |2|b in the range of masses between 1.5 and 1.7 GeV is mainly due to the 
decay D° —> K~n + n °, where the 7T° is not measured. A fit of the M(Ktt) spectrum in the mass 
range 1.7 < M(Ktt) < 2.5 GeV to an exponential curve plus a Gaussian gives a value for the 
D° mass of 1858 ± 3MeV, slightly below the PDG|33| value of 1864.6 ± 0.5MeV. The width 
of the signal is 19 ± 3 MeV. The fit result is shown by the solid curve in Fig. |2b. The signal 
region in AM extends from 143 to 148 MeV. The background under the signal is estimated in 
two independent ways: 


• counting the number of combinations which use pairs of tracks with the same charge 
for the D° pre-candidate and fulfill the above requirements (wrong charge combinations, 
shown in Fig. ^]a as the solid histogram); 

• using a control region 2.0 < M(Ktt) <2.5 GeV instead of the signal region 1.80 < 
M(Ktt) <1.92 GeV where the events from the control region have been normalized to 
the number of events with 155 < AM < 180 MeV. 


The number of background candidates is taken as the weighted mean of the two estimates. The 
number of reconstructed D* candidates is taken to be the number of entries in the signal region 
minus the number of background candidates. For 5 < Q 2 < 100 GeV 2 , y < 0.7 in the restricted 
kinematic region 1.3 < pt(D*) < 9.0 GeV and \r](D*)\ < 1.5 the result is 122 ± 17 D* mesons 
above a background of 95 ± 8 . The AM and D° mean value and width of the signals predicted by 
the NC DIS MC sample agree well with the data. The same is true for the background shape 
and the signal to background ratio. For the measurement of the p^(D*), r](D *), Q 2 and W 
dependence, the data sample in each variable has been divided in three bins. The number of 
D* obtained in each bin is determined following the same procedure described above. 


5.2 D* reconstruction efficiency 


The D* reconstruction efficiency is determined using MC simulations. DIS events with cc pro¬ 


duction by photon-gluon fusion were generated using two different MC models: AROMA 2.1 134 
and HERWIG 5.8 |35|. AROMA is a MC model for heavy flavour production. It is based on the 
following ingredients: (i) the complete matrix elements in LO for the PGF process 7 * g —> cc 
(taking into account the mass of the charm quark and the full electroweak structure of the 
interaction), ( ii ) gluon emission from the cc system in a parton shower approach, (Hi) initial 
state parton showers, ( iv ) hadronisation with the Lund string model [^Dj as implemented in 
JETSET 7.4 [|30|] and heavy flavour decay. In addition bremsstrahlung emission of gluons (some 
of which may split perturbatively into two gluons or qq pairs) is also included. 


HERWIG is a general purpose QCD MC event generator for high energy hadronic processes. 
Here it is used for simulating 7 * g —> cc in LO. I 11 addition, leading-logarithm parton show¬ 
ers were included in the simulation. Fragmentation into hadrons is modeled with a cluster 
fragmentation model which takes into account the charm quark mass. 
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The parameters of the MC programs were set to their default values. In particular, the charm 
quark mass (■ m c ) was set to 1.35 GeV for AROMA and 1.8 GeV for HERWIG. For the parton 
densities of the proton the parameterisations of MRSD—' |32j and MRS A were considered. 


Those events containing at least one charged D*, decaying into D° 7 r + with subsequent decay 
D° —> K~tt + , were processed through the standard ZEUS detector and trigger simulation 
programs and through the event reconstruction package. Approximately 5000 MC events from 
the samples generated and processed in this way, passed all the selection criteria. Both MC 
reproduced equally well the shapes of the uncorrected data distributions relevant for this study. 
These events were used to determine the D* reconstruction efficiency as a function of p^(D*) and 
7](D*). For each D* candidate, the efficiency in a given (pt(D*), r](D*)) bin is the ratio of the 
number of reconstructed D* in the bin to the number of generated D* in the bin. Both 
AROMA and HERWIG MC samples give compatible results and were combined for the final 
efficiency determination. The r](D*) resolution is less than 0.01 in units of pseudorapidity. 
The pt(D*) resolution is about 2% of the bin width chosen for the figures; none of these 
show systematic shifts. The D* reconstruction efficiency in the kinematic range considered 
varies between 20% for low p^(D*) and 70% for p^(D*) >5.8 GeV. The efficiency varies in r] 
from 30% near | r] \ = 1.5 to 50% for D* moving transversely to the beam direction (jj = 0). 
The average D* reconstruction efficiency is about 38%. As described in section 4.3, the event 
selection efficiency is determined using the HERACLES MC. The convolution of the event 
selection efficiency and the D* reconstruction efficiency gives an overall detection efficiency of 
approximately 30% for D* decaying via (1). 


6 Results 


6.1 D* fractional momentum distribution 


To investigate the charm production mechanism, the distribution of the fractional momentum of 
the D* in the system [xd* = lt % ) is studied. The PGF process produces a cc pair that, in 
the 7 *p system, recoils against the proton remnant. In contrast, in the flavor-excitation process, 
a single c quark from the proton sea is scattered off the proton, flying, in the system, in a 
direction opposite to the proton remnant. Since, in general, the D* carries a large fraction of 


the momentum of the parent charm quark [[k|], clear differences are expected between the Xd* 
distributions from the two production mechanisms. The distribution is expected to be centered 
at xd * < 0.5 for PGF and peaked at high xd* values for direct production. Figure [l]c shows 
the normalized xd* differential distribution measured in the data and the AROMA prediction 
(based on PGF). For comparison we also show the prediction by LEPT06.1 with only direct 
production from charm. Both MC use the proton density of MRSD—' |32|. The shape of our 


data distribution is compatible with PGF, in accord with the HI result |14 . 


6.2 Cross sections in a restricted pt(D*) and rj(D*) kinematic region 

Differential cross sections as well as cross sections integrated over the kinematic region 1.3 < 
Pt(D*) < 9.0 GeV and \rj(D*)\ < 1.5 are presented in this section. The e + p —> e + D* ± X 
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differential cross sections are corrected for the efficiencies of the selection criteria as well as for 
the branching ratios B(D* + —■> D°irf) x B(D° —► K~n + ) = 0.0262 ± 0.0010 

Fig-0 shows the differential e + p —> e + D* ± X cross section as a function of a) pr(-D*), b) r](D*), 
c) W and d) Q 2 in the kinematic region defined above. The inner error bars are the statistical 
errors and the outer ones show the statistical and systematic errors added in quadrature. The 
data points in the p^(D*) and Q 2 distributions are shown at the positions of the average values 
of an exponential and polynomial fit, respectively, for a given bin. 

The D * ± differential cross section da / dpr 2 (D*) exhibits an exponential falloff in px{D*) and 
da/dr](D*) is approximately flat. The shape of the W dependence of the cross section predicted 
for photon-gluon fusion at low W is determined by the p^{D*) > 1.3 GeV and the \r](D*)\ < 1.5 
cuts. The falloff at large W is mainly due to the virtual photon flux. The D * ± cross section 
drops steeply as Q 2 increases. 



The systematic uncertainties arise from several effects. The main uncertainties coming from 
DIS event selection and D* reconstruction are: the transverse momentum cut of the soft pion 
(7%), the range used for the K tt invariant mass (5%) and the positron energy cut (5%). The 
differences in acceptances evaluated from different MG samples correspond to an uncertainty 
in the integrated e + p —» e + D* ± X cross section at the 7% level. The QED radiative effects 
contribute 3 % to the systematic error. The systematic error attributed to the branching ratios 
is 4%. The normalisation uncertainty due to the determination of the luminosity and trigger 
efficiency is 2%. No significant contribution to the systematic error from the primary vertex 
requirement was found. Adding the contributions in quadrature, a total systematic error of 15% 
is obtained in the determination of the integrated cross section. For the differential distributions 
the systematic errors are estimated bin by bin. 


The results of the NLO analytic calculation of |D] using the GRV|37| NLO gluon density, 
xg(x, Q 2 ), of the proton is shown as a band in Fig. where the upper (lower) limit corresponds 
to a charm quark mass of 1.35 (1.7) GeV. The calculation is based exclusively on the PGF 
process and no charm produced from the proton sea is considered. The GRV NLO parton 
density gives a reasonable description of the proton structure function F 2 measured at HERA 
[£i], |2l[ and treats the charm quarks as massive particles, consistent with the ansatz of |lj. 
Charm quarks were hadronized according to the Peterson fragmentation function with 
e c = 0.035 ± 0.009 |I9| and the probability for a charm quark to fragment into D* , P(c —> D*), 
was assumed to be 0.26 ±0.02 |39|]. Both the renormalisation and the factorisation scales were 


chosen to be g = \JQ 2 + 4 m 2 . In this calculation the largest uncertainty is related to that of 
the charm quark mass. The predicted differential cross section for D*± production as a function 
of a) Pt(D*) and b) rj(D*) is shown in Fig. |3| outside the restricted kinematic region selected for 
the data. The NLO analytic calculation reproduces well the shapes of the p T (D*),W and Q 2 
distributions and is consistent with the r](D*) dependence in the restricted kinematic region. 


Below, unless explicitly stated otherwise, the LO and NLO analytic calculations are those from 
[|TJ and use respectively the LO and NLO gluon density parameterisations of GRV j37J and a 
charm quark mass m c = 1.5 GeV. 


In Fig. the data are also compared with the LO analytic calculation, shown as histograms. 
The LO calculation performed with the MC programs used in the acceptance calculation and 
described in section 5.2 agrees well with the LO analytic calculation. Note that the MC 













calculation includes parton showers in addition to the LO matrix elements. The LO MC and 
LO analytic calculations describe the shapes of the p^(D*) , Q 2 and W dependences equally 
well. 

Integrated over the region 5 < Q 2 < 100 GeV 2 and y < 0.7, the e + p —> e + D* ± X cross section 
for D * ± in the restricted kinematic range 1.3 < pt{D*) < 9.0 GeV and \r](D*)\ < 1.5 is 

ap T ,v( e+ P e + D*± X) = 5.3 ± 1.0 ± 0.8 nb 

where the first error is statistical and the second error is systematic. The e + p —» e + D* ± X 
cross section represents about 4 % of the DIS cross section (calculated from structure functions 
determined by fits to global data |32|) in the same Q 2 and y ranges. Tabled compares the 
measured cross section with the LO and NLO predictions. The integrated cross sections pre¬ 
dicted in NLO is in reasonable agreement with the data. The LO MC and analytic calculation 
agree with the data, although their sensitivity to the parameters involved in the program (m c , 
xg(x, Q 2 ), e c and /i) is higher than in the case of NLO calculations. 


6.3 Integrated charm cross section and F£ c 

The DIS inclusive cross section for charm production, e + p —> e + ccA", <j cc , can be expressed in 
terms of F 2 CC using: 


d 2 a c 


2i ra 2 


dxdQ 2 Q 4 x 


[(1 + (1 - y) 2 ) Ff(x, Q 2 ) - y 2 Ff{x, Q 2 )} 


( 2 ) 


In the standard model, the contribution from Z-boson exchange is expected to be small in the 
( Q 2 ,x ) range of the present analysis and therefore the F| c contribution is neglected. The F£ c 
contribution has been estimated m to be smaller than 2% and, therefore, no correction has 
been applied. 


In order to estimate cr cc and to evaluate F| c as a function of Q 2 and y, the measurements must 
be extrapolated to the full p T (D*) and g(D* 
cross section using: 


cr(e p 


ccX) = 


) range. a cc is obtained from the integrated D 
1 a(e + p —> e + D* X) 


(3) 


2 P(c —> D* 

According to the extrapolation with the NLO calculation, about 50 % (at Q 2 = 45 GeV 2 ) to 
65 % (at Q 2 = 7 GeV 2 ) of the D* production is found to be outside the restricted kinematic 
region. Using the LO analytic calculation for the extrapolation results in similar cross sections. 

For the determination of the integrated cross section the data are divided into two Q 2 in¬ 
tervals, namely 5 < Q 2 < 10 GeV 2 and 10 < Q 2 < 100 GeV 2 . The extrapolated accep¬ 
tance is calculated using the NLO prediction of GRV[[D]] with m c = 1.5GeV,// = \JQ 2 + 4m 2 , 
P(c -> D* + ) = 0.26 ± 0.02 and e c 

tion, cr (e + p —> e + ccA"), for y < u.i, ia _i_ _i_ ±.o _ x 2 

12.5 ± 3.1 ± 1.8 jrilfnb for 10 < Q 2 < 100 GeV 2 . The first error gives the statistical, the see 


= 0.035 ±0.009 |39[]. The integrated charm cross sec- 


0.7, is 13.5 ± 5.2 ± 1.81];'® nb for 5< Q 2 <10 GeV 2 , and 


ond the experimental systematic uncertainty and the third the model dependent uncertainty. 
The latter is studied by varying the relevant parameters of the model used for extrapolation, 
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namely m c (from 1.35 to 1.7 GeV), /j (from 2 m c to 2 y^Q 2 + m 2 ), e c (from 0.025 to 0.045), and 
xg(x, Q 2 ) (GRV [|37H , MRSG [[HU , CTEQ3 |42|1 , MRSA'^BH) [] and also includes the error in the 
probability for c —> D* . An uncertainty between —12% and +9% is obtained for both Q 2 
ranges. 


The integrated charm cross sections as well as the corresponding predictions from the LO 
and NLO analytic calculations are also listed in Table [I]. The NLO predictions are about one 
standard deviation smaller than the measured value in the first Q 2 interval (from 5 to 10 GeV 2 ) 
but in agreement in the second Q 2 interval (from 10 to 100 GeV 2 ). The charm cross section 
obtained for 10 < Q 2 < 100 GeV 2 is in reasonable agreement with the recent HI measurement 
in the same Q 2 region. 



a PT:V [nb] 

(T 1 [nb] 

og [nb] 

Data 

5.3 ± 1.0 ± 0.8 

13.5 ±5.2 ±1.8 7^ 

12.5 ± 3.1 ± 1.8 ij;? 

AROMA 

4.57 

12.6 

14.2 

LO Analytic calculation 

4.79 

11.0 

12.4 

NLO Analytic calculation 

4.15 

9.4 

11.1 


Table 1: The measured e + p —> e + D* ± X cross section for 5 < Q 2 < 100 GeV 2 , y < 0.7 in the 
restricted kinematic region 1.3 < pt(D*) < 9.0 GeV and \y(D*)\ < 1.5 ( ^p T ,rj ); the extrapolated 
e + p —> e + ccX cross section for y < 0.7 and for 5 < Q 2 < 10 GeV 2 (ay) and 10 < Q 2 < 100 GeV 2 ( 02 ). 
The measured cross sections are compared with the NLO analytic calculation for the GRV NLO gluon 
density as well as with the LO MC and LO analytic calculations. The statistical uncertainties of the 
NLO analytic calculation have about 0.01 the size of the uncertainties of the data points. 


Fig. [I] shows the resulting F 2 CC values as a function of x for the different Q 2 and y bins as defined 
in Fig. 0b. The HI measurements |14[ are also shown; they were taken in the same Q 2 bins. 
The two sets of data are in good agreement. Also included in Fig. [J are the data from the 
EMC fixed target experiment H which were measured at x values between 0.02 and 0.3. The 
charm contribution, i% c , to the proton structure function is seen to rise by about one order of 
magnitude from the high x region covered by the fixed target experiments to the low x region 
measured by HI and ZEUS. The NLO analytic calculation for F| c is shown in Fig. [4| as a band 
(the upper (lower) limit corresponds to a charm quark mass of 1.35 (1.7) GeV). In the region of 
our measurements, the use of CTEQ4F3 ]44| or MRRS ||5]] gluon densities give results which 


are within 5 % of those using GRV. In the present Q 2 range, F 2 C scales roughly with the input 
gluon density in NLO perturbative QCD. The measured rise of F 2 CC from the high x to the low 
x region is reasonably described by NLO perturbative QCD. 


Comparing with our measurement of F 2 (x,Q 2 ) El , we observe that the ratio F| c /F 2 is about 
25% for the entire (Q 2 ,x) range of the present analysis. 


3 Note that while GRV (and also CTEQ4F3 and MRRS, see below) follow the charm quark treatment of [Q 
based on PGF, the CTEQ3 and MRS distributions treat charmed quarks as massless quarks with a parametrized 
density contributing above a certain threshold in Q 2 and therefore are not completely consistent with the 
program of |I[. 
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7 Summary and conclusions 


We have measured the D* differential and integrated e + p —> e + D* ± X cross sections for deep 
inelastic scattering at a/s = 300 GeV with 5 < Q 2 <100 GeV 2 and y < 0.7 in the restricted 
kinematic region of 1.3 < p^(D*) <9.0GeV and \r](D*)\ < 1.5. The integrated D cross 
section is measured to be 5.3 ± 1.0 (stat.) ± 0.8 (syst.) nb. 

The shape of the D* fractional momentum distribution in the rest system, Xd*, shows that 
the PGF mechanism prediction agrees well with the data for DIS charm production in this 
kinematic range. 

A QCD analytic calculation in NLO with the NLO GRV gluon density reproduces the shapes 
of the pt(D*) W and Q 2 distributions and is consistent with the p(D*) dependence in the 
restricted kinematic region. The predicted cross sections are in reasonable agreement with the 
data. We have used QCD calculations to extrapolate the D cross section measured in the 
restricted pT(D*),r](D*) region to the full region and estimated the integrated charm cross 
section and the charm contribution F 2 C to the proton structure function F 2 . When compared 
to the fixed target measurements (performed at large x ) F 2 C is found to rise as x decreases. The 
rise is described by NLO QCD calculations when using a gluon density consistent with that 
extracted from the scaling violations in the proton structure function F 2 measured at HERA. 
Such a gluon density distribution is also compatible with our previous measurements of D* in 
photoproduction. 
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X D* 

Figure 1: (a) LO diagram for photon-gluon fusion (PGF). (b) The (x,Q 2 ) plane with the Q 2 , 

y region and bins chosen for the Ff c analysis. Large dots correspond to the D* candidates; 
small dots correspond to a subsample of DIS events. The second constant y line corresponds to 
y = 0.2. (c) Normalized e + p —> e + D* ± X cross section for 5 < Q 2 < 100 GeV 2 , y < 0.7 in the 
restricted kinematic region 1.3 < p^(D*) < 9.0 GeV and \rj(D*)\ < 1.5 as a function of xr>*. 
The inner error bars show the statistical errors and the outer ones the statistical and systematic 
errors added in quadrature. The horizontal bars represent the bin widths. The prediction for 
PGF as calculated with AROMA (solid histogram) and the charm sea contribution as calculated 
with LEPTO 6.1 selecting QPM events (dashed histogram) are also shown (see text). 



14 
















AM [GeV] 



M(Kk) [GeV] 

Figure 2: (a) AM mass distribution for Kir combinations in the M(Kn) signal re¬ 

gion (1.80 < M(Kir) < 1.92GeV), full dots, and for the wrong charge combinations, 
solid histogram, (b) M(Kn) mass distribution for the Kirir s combinations in the AM = 
M(Kirir s ) — M(Kir) signal region (If 3 < A M < lf8MeV), full dots, and for the wrong 
charge combinations, solid histogram. The solid lines in both figures show the result from the 
fits (see text for details). 
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Figure 3: Differential e + p —> e + D* ± X cross sections for 5 < Q 2 < lOOGeV 2 , y < 0.7 in the 
restricted kinematic region 1.3 < p^(D*) < 9.0 GeV and \r)(D*)\ < 1.5 as a function ofp^fD*) 
(a), r)(D*) (b), W (c) and Q 2 (d). The inner error bars show the statistical errors and the outer 
ones correspond to the statistical and systematic errors added in quadrature. The horizontal 
bars represent the bin widths. The NLO QCD prediction for different charm quark masses is 
shown by the band (see text). The LO prediction for the GRV(LO) gluon density is shown by 
the histogram (see text). The predicted cross sections in (a) and (b) are shown without the 
p T (D*) and rj(D*) cuts respectively. 
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Figure 4: The charm contribution, Ff c , to the proton structure function F 2 as derived from 
the inclusive D*(ZEUS and HI) mid D°(H1) production compared with the NLO QCD predic¬ 
tions based on the GRV parton distribution using different charm quark masses for Q 2 =7, 12, 
25 and f5 GeV 2 The upper (lower) limit of the band corresponds to a charm quark mass of 
1.35 (1.7) GeV (see text). The results from the EMC collaboration are shown as crosses. For 
the ZEUS data, the inner error bars show the statistical errors and the outer ones correspond 
to the statistical and systematic errors added in quadrature. The error bars from HI show the 
statistical and systematic errors added in quadrature. The error bars for EMC are within the 
symbol. 
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